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ABSTRACT 

Ambient interstellar material may become entrained in outflows from massive 
stars as a result of shear flow instabilities. We study the linear theory of the 
Kelvin - Helmholtz instability, the simplest example of shear flow instability, 
in a partially ionized medium. We model the interaction as a two fluid system 
(charged and neutral) in a planar geometry. Our principal result is that for much 
of the relevant parameter space, neutrals and ions are sufficiently decoupled that 
the neutrals are unstable while the ions are held in place by the magnetic field. 
Thus, we predict that there should be a detectably narrower line profile in ionized 
species tracing the outfiow compared with neutral species since ionized species 
are not participating in the turbulent interface with the ambient ISM. Since the 
magnetic field is frozen to the plasma, it is not tangled by the turbulence in the 
boundary layer. 

1. Introduction 

In the process of forming, massive (> 10 Mq) stars produce massive bipolar outfiows. 
In contrast to outfiows from low-mass proto-stars, outfiows from massive stars frequently 
have more mass than their presumed driving proto-star (Churchwell 1997). The canonical 
explanation for the origin of the mass is that a less massive jet originating from either the 
accretion disk or protostar sweeps up and entrains ambient interstellar material. This model 
can potentially explain both the large masses and poor collimation of massive outfiows. It 
may also explain how low-mass outfiows become less coUimated as they age (Bachiller and 
Tafalla, 1998). Churchwell (1997) argued that mechanical entrainment, through e.g. bow- 
shocks, could not simultaneously explain both the large mass and young age of massive 
outfiows. 
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Shear flow, or Kelvin- Helmholtz, instabilities at the interface between the jet and the 
ambient material are an alternative possible mechanism for entrainment. The instability 
amplifies ripples at the interface, culminating in the development of a turbulent boundary 
layer which can transfer momentum from the jet to the ambient medium. A magnetic 
field parallel to the fiow tends to suppress the instability because magnetic tension opposes 
rippling. The result is that instability is only present if the velocity shear is greater than the 
Alfven velocity (Chandrasekhar, 1961). 

The nonlinear development of the Kelvin-Hclmholtz instability is strongly infiuenced by 
the magnetic field (Malagoli et al., 1996 & Jones ct al., 1997). The numerical simulations 
performed by these authors show that the vortices created by the instability wind up the 
field until it reaches the Ohmic scale. Windup transfers energy to small scales, while Ohmic 
dissipation heats the layer. These effects are absent in the field-free case. 

In the following study, we investigate the Kelvin- Helmholtz instability in the linear, 
partially ionized regime to determine its possible contribution to entrainment in massive 
bipolar outflows. This problem cannot be analyzed in the hydrodynamic regime because 
the neutral Alfven velocity (v^ = ^^^^ ) is sometimes larger than the flow speed, indicating 
that the magnetic held can be strong enough to affect the dynamics. On the other hand, 
the pure magnetohydrodynamic (MHD) regime is inappropriate because at sufficiently short 
wavelengths, Alfven waves propagate nearly independently of the neutrals. 

In §2 we describe a model that we believe is appropriate for massive outflows driven by 
massive protostars. In §3 we discuss the numerically-derived solutions to the dimensionless 
dispersion equation and equivalent analytical solutions to simplified versions of the dispersion 
equation. We demonstrate that there is a new instability, unique to the partially ionized 
regime, which results from ion-neutral slip. In §4, we discuss the observational implications 
of our results, concentrating on the physical environment of the interaction zone between a 
bipolar outfiow and the ambient interstellar material (ISM). §5 is a critical summary of the 
paper. 



2. Model Description 

We model the interaction zone between a high speed jet and the ambient ISM by a slab 
geometry with a fiow along the x-axis. The interface between the two fiuids is the x — y 
plane. The density is everywhere uniform, and the velocity is uniform except at 2; 0, where 
it has a step. We assume the flow velocity U decreases with increasing z. The magnetic held 
is uniform and aligned with the flow, which maximizes its stabilizing effect in the x-direction. 
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In this model, the equihbrium flow has zero vorticity except on the plane z — 0, where 
the vorticity is infinite. For our choice of coordinates and velocity profile, the total vorticity 
integrated across the sheet is —2Uy. 

The mechanism of the Kelvin-Helmholtz instability is explained in Batchelor (1967). 
Suppose the interface is perturbed by a small amplitude ripple, as shown, in the center of 
mass frame, in Figure 1. Denote the vertical displacement of the vortex sheet by 77. In 
general, the vorticity perturbation is out of phase with 1]. If it can be arranged that the 
vorticity has the distribution shown by the dashed line, then the equilibrium flow sweeps 
negative vorticity toward point P on both sides of the interface. The effect is as if an extra 
eddy, rotating counterclockwise as sketched in the Figure, were added to the flow. The 
centrifugal force associated with the eddy pushes the interface between P and Q outward, 
resulting in the growth of 77 with time. If the vorticity perturbation were shifted in phase 
by TT with respect to the dashed line, 77 would shrink. As explained here, the instability 
mechanism requires advection of vorticity by the bulk flow in opposite directions above and 
below the interface. It is always possible to choose a reference frame in which this is the 
case, provided that the velocity profile has a jump. 

A magnetic field parallel to the fiow exerts a tension force when bent. The tension 
force is stabilizing, and the criterion that it dominates hydrodynamic forces is roughly that 
the magnetic energy exceeds the kinetic energy of relative motion, or, equivalently, that the 
Alfven speed for the bulk medium exceed the fiow speed. 

The basic mechanism of instability is unaffected by partial ionization. In the weakly 
ionized medium of interest here, the plasma component tends to be stable, because its low 
density implies low kinetic energy. The neutral component, on the other hand, is intrinsically 
unstable because it is not directly acted upon by magnetic forces. Friction between the ions 
and neutrals allows magnetic forces to act on the neutrals. As we will see, however, no 
matter how large the magnetic field is, it cannot stabilize the neutrals completely. 

In the present discussion, we follow the analysis of Chandrasekhar (1961) but consider 
two ffuids, charged and neutral, coupled by elastic collisions^. We assume the medium is 
incompressible (V • Vj „ = 0). Our problem is described by the neutral and ion momentum 
equations and the magnetic induction equation: 



We ignore coupling by ionization and recombination, since these processes are slow in comparison to 
elastic processes. 
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PidtVi^ PiVi-Vwi = -VPj -p^i/niK - v„) + ^(V X B) X B (2) 
dtB = Vx(viXB), (3) 

where 

l^ni = IPi (4) 
is the neutral-ion coUision frequency, and 

7= 5 

is the coUision rate coefficient per unit mass. 

The system consisting of two constant density, constant pressure regions of uniform 
flow ahgned with a uniform magnetic field is clearly a steady state solution of eqns. (1 - 
3). We now consider small perturbations of the fiow. We assume these perturbations are 
independent of y. According to Squires' Theorem (Drazin & Reid 1981) and its extension 
to MHD (Hughes & Tobias 2001), in the absence of viscosity and resistivity these are the 
fastest growing perturbations. Since the magnetic field enters the problem only through 
the tension force, which is independent of the structure in ion-neutral friction should not 
affect this result. We expand the momentum and induction equations to first order in a 
small parameter e, writing 

Vx = u — U -\- eu (6) 
Vy = w — ew (7) 
Bx = B + ebx (8) 
B, = eb,. (9) 

We assume all perturbed quantities q take the form 
where qq is a constant amplitude. Linearizing in e gives 

pn{dt + Ud^)Un = -d^Pn + pn^niiUi - lln) (U) 
Pn{dt + Ud^)Wn = -dxPn + Pnl^ni{Wi - Wn) (12) 
pi{dt + Udx)Ui = -dxPi + pnVni{Un - Ui) (13) 

Pi{dt + Udx)Wi = -dxPi + pn'^niiWn - Wi) + —(dxb^ - d^bx) (14) 

{dt + Udx)bx = BdxUi (15) 
(dt + Udx)b, = BdxWi, (16) 



- 5 - 



which correspond to the equations (180-183) in Ch. 12 of Chandrasekhar (p. 508-9). Ehm- 
inating B in equation (14) using equations (15) and (16), constructing the z component of 
the curl of the momentum equations, and applying V • Vj_„ = 0, we obtain 

Pi{cU + UK) - _ , , ) (dl - kl)Wi = %Pnl^ni{dl - kl){Wi - Wn) (17) 



Pn{uj + Uk^){d^-kl)wn = ipnVrii{dl-kl){wn-Wi) (18) 

These relations hold in each region, but not at the interface, where the velocity is not analytic. 
Thus, to derive the dispersion relation, we must integrate each equation over an infinitesimal 
region surrounding z — Q and use standard Gaussian pillbox arguments. Denoting the jumps 
in quantities by A{}, we get 

^{pi{u + Uk:,)d,Wi} = i^{pnVnid,{Wi-Wn)} + ^-^^{ ^""^^ } (19) 
/^{Pn{uJ + Ukx)d^Wn} = i/^{pnynidz{Wn-Wi)}. (20) 

We write the ion and neutral velocity amplitudes in the form 

Wii = Aicu + Uik)e''^ z<0 (21) 

= A{u; + U2k)e-''' z>0 (22) 

Wni = B{u; + Uik)e^^ z<Q (23) 

Wni = B{uj + U2k)e-^' z>0, (24) 

where the subscripts 1 and 2 refer to the jet and ambient medium, respectively. The forms 
of Wi^n are chosen such that the ion and neutral displacements are continuous across the 
interface, as must hold on physical grounds. Substituting these forms into equations (19) 
and (20), we solve for B, the neutral velocity fluctuations, in terms of A, the ion velocity 
fluctuations. The result is 

^ _ _^ PniZ/i(cc; + Ulkx) + Pn2l^2{uJ + U2k^) ^ 

p„i(a; + Uikx){ixi + Uikx - ivi) p„2(<^ + U2kx){uj + U2kx - ii'2) 
We then use equations (21) and (25) in equation (19) to obtain the dispersion relation 

Pii{u + Uikf + Pi2{uJ + U2kf = {pniM^ + Uik) + Pn2l^l{uJ + U2k)) X (26) 

{pnlT^li'^ + Uik) + pn2T^l{u + U2k)) 

X — I H~ 

Pni{oj Uik){uj Uik — iui) + Pn2{uJ + U2k){uj U2k — ii'2) 

B^e 

f 



27r 

where i/ni has been abbreviated to u and k^ abbreviated to k. 
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We now assume that the jet and ambient medium have the same density, work in the 
center of momentum frame so that 

(27) 





= -U2 = 


u 


Pnl 


= Pn2 = 


Pn 


Pil 


= Pi2 = 


Pi, 



and nondimensionahze the problem as follows: 



m = ^ (28) 

Pi 
UJ 

X = — 

V 

h = — 

V 

^/mBk kvAi 
a — 



v^fWfha y 

where the ion Alfven speed vai is defined as The variables ,t, h and a represent the 

perturbation frequency, flow speed and ion Alfven speed, respectively. 

Substituting equations (27) and (28) into equation (26) we obtain 

-i{m^\)x^ ^{2h^ - a^)x^ -%{mh^ - ^h^)x^h^{h^ - a^) = 0. (29) 

In all cases of interest, m 3> 1 and a ^ h, which allows us to drop a few of the terms 
in equation (29). This leads to the somewhat more compact expression 

- imx^ - o?x^ - i{mh^ - a})x - a^h"^ = 0. (30) 

Equation (30) is the basis of our subsequent stability analysis. We will be looking for solutions 
for X which have a negative imaginary component, an indication of a growing instability. 



3. Analysis of the Dispersion Relation 

Equations (29) and (30) reduce, in special cases, to previously known results. These 
arc discussed in §3.1. In §3.2 wc present numerical and analytical solutions for the fastest 
growing modes. In §3.3 we discuss the physical nature of the instabilities. 
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3.1. Recovery of previously known results 

If h = 0, equation (30) reverts to the dispersion relation for Alfven waves in a partially 
ionized medium (Kulsrud & Pearce 1969) 

X {x^ - a^) ^ [x- - ^) = 0. (31) 

At sufficiently small wavenumber {k < a < the wave speed is close to the neutral 

Alfven velocity vau because the neutral-ion collision timescale is less than the wave period 
and the ion and neutral species are closely coupled. We can derive this result from equation 
(31) by treating the second factor in parentheses as dominant. This leads to 

a , , 

x^±^ + i—, 32 



which is a weakly damped Alfven wave propagating at speed vau- 

At sufficiently large wavenumber {k > a > m/2), the wave period is shorter 

than the ion-neutral collision timescale. Thus, Alfvenic disturbances are able to propagate 
through the ions without disturbing the neutrals and propagate at the ion Alfven velocity 
VAi- We can derive this result from equation (31) by treating the first factor in parentheses 
as dominant. This leads to 

xp=i±a + i—, (33) 

which is a damped wave propagating at the ion Alfven velocity. In the intermediate wavenum- 
ber regime (2^Jrn < a < m/2), waves cannot propagate because the ions are poorly coupled 
to neutrals and thus feel a strong drag. 

The shear instability in a fully-neutral medium (or ionized but magnetic field free) is 
present at any velocity and has growth rate kU (= h in our notation). In order to recover 
the hydrodynamic instability we must use equation (29) with a = 0. The dispersion relation 
can then be written in the form 

(x^ + h'') [x^ -iim + l)x + h''] = 0. (34) 

The hydrodynamic instability is the root x = —ih. 

On the other hand, a fully ionized medium with a magnetic field parallel to the shear 
fiow is stabilized by magnetic tension if U < va- In order to recover the MHD instability, 
we set m = in equation (29) and write the dispersion relation in the form 

{x^ + h^- a^) {x^ + h^- ix) = 0. (35) 
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Note that the growth rates of both hydrodynamic and MHD instabihties increase hnearly 
with wavenumber. 



3.2. The most unstable roots 

With §3.1 as background, we now turn to the shear instabihty in the partially ionized 
regime. 

If there were no ion-neutral friction, the neutrals would always be unstable. The plasma 
would be stable iiU < vai {h < a), which we assume to be the case. Since the growth rate of 
the hydrodynamic instability increases with there is a wavenumber above which friction 
with the ions cannot stabilize the perturbation. Therefore, we expect short wavelength 
perturbations to grow at the hydrodynamic rate in the neutral fluid, but to leave the plasma 
and magnetic field in place. 

At longer wavelengths, there is time for ion-neutral friction to act. Hva <U {a < hy/m), 
the magnetic field cannot be expected to stabilize the medium no matter how well the ions 
and neutrals are coupled. On the other hand, iiU < va, we might expect that at sufficiently 
long wavelengths the coupling is good enough to completely stabilize the system. We will 
see, however, that this is never the case. 

We carried out a parameter study of the roots of equation (30) as functions of a and 
h using the Mathcmatica software package (Wolfram 1999). Surprisingly, the roots are 
closely approximated by the analytical expressions given in Table 1 for the specified range 
of dimensionless parameters. In the results presented here, we focus on two sets of physical 
environments which we believe arc illustrative of the range of conditions encountered in 
molecular outflows. The density of the ambient medium near massive protostars is a strong 
function of distance from the protostar. within radii of 100-1000 AU densities are ^10^ 
cm~^, further out where outflows can be easily resolved and structure studied densities of 
10^ cm~^ are typical. In both examples, we take the molecular hydrogen density n{H2) to 
be 10^ cm~^, the relative flow speed 2U to be 20 km s~^, and the ratio m of neutral to 
ionized mass density to be 10^ (corresponding to an ionization fraction ~ 10~^ for an ion 
to neutral particle mass ratio ~ 10). The collision frequency v is 1.5 x 10^^^ s~^ (we use 
the rate coefficient of Draine, Robcrgc, & Dalgarno (1983); {av) = 1.5 x 10^^ cm^ s^^). 
With the these values, the wavelength at which the hydrodynamic frequency kU matches 
the coUison frequency {h = 1) is about 1.3 pc. The only difference between the two cases is 
in the magnetic fieldstrength, which we take to be 1 mG (strong field case) and 1 /iG (weak 
field case). With these parameters, a/h — 4400 and 4.4 in the strong and weak field cases. 
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respectively. 

Since only growing modes will be dynamically important in the interface between a fast 
moving jet and the ambient ISM, we restrict our discussion to two solutions, —ih [and the 
variant —i{^ — h)] and —ih^. Confirmation that these expressions closely approximate the 
numerically derived roots is shown in Figure 2. The numerically determined roots to the 
dispersion equation arc shown as a function of wavelength; overplotted are the analytical ap- 
proximations to the two growing modes. The wavelengths 27r/A; and growth times {Im{u}))~^ 
are plotted in physical units for the two cases specified above. 

The analytical solutions can be derived from the dispersion relation by dropping small 
terms. The dominant terms used to derive each solution are given in Table 1. For example, 
the first root listed, ~ih, is the standard hydrodynamic instability derived from equation 
(34). From the dispersion relation, we can see that the terms oc (— a^x^, ia^x, —a^h"^) 
dominant in the high magnetic field regime. Of the terms, the terms oc /i^, dominate 
over the term oc x because h> 1. In contrast, when h <1, the terms cc x, h'^ dominate, and 
the growing mode is —ih^. 

This second growing mode, —ih"^, is a new mode that arises only in a weakly ionized 
medium. It originates when /i < 1, or equivalently kU < v. The neutrals feel a constant drag 
force, independent of position, caused by collisions with the stationary ions. The instabihty 
has an origin similar to the standard hydrodynamic instability explained in §2. The principal 
difference is that the neutral-ion collision time is shorter than the wave period. Thus, neutral 
material collides with ions (held stationary because of the magnetic field tension) frequently 
during each perturbation period, which slows the growth of small perturbations. This slowing 
of the growth is clearly visible in Figure 2 (bottom panel) where the growth time is longer for 
long wavelength (and long period) perturbations. The origin of this mode can be made more 
clear by assuming the ions are held strictly in place by the magnetic field and re-deriving 
the roots (see Appendix A). The break between the two growing modes occurs when h = 1 
or, equivalently, ^ (the wave period) = ^ (the neutral- ion collision timescale). See Figure 3 
for physical parameters of this transition. Since this mode grows more slowly than the —ih 
mode, it is important only in the situation U < va (or ^/mh < a), when the MHD mode is 
stable. 

When ^/rnh > a and h < 1, —ih-\Jl — a'^/mh'^ is a better approximation to the growth 
rate. This is the analog of eqn. (35) for a fully ionized medium. It is based on the limit of 
infinitely strong ion-neutral coupling. 
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3.3. Physical nature of the unstable modes 



One of the main results of §3.2 is that the plasma and the neutrals can follow very 
different dynamics as the instability grows. Although the interface between the two regions 
of flow is the same for the plasma as for the neutrals in the equilibrium system, the dis- 
placed interfaces for the plasma and neutral species can differ, as can their velocities. The 
negligibly small inertia of the ions renders any differences irrelevant to the basic mechanism 
of instability, but they are important for the nonlinear outcome of the instability, and for its 
spectroscopic signatures. 

On physical grounds, we expect large differences between Vi and f„ when ion - neutral 
couphng is weak {h > 1) and magnetic forces are strong. Whether these differences persist 
at small h depends on the relative values of aj^fm and h. If their ratio is less than unity 
the field is too weak to suppress instability even in the limit of infinitely strong ion-neutral 
coupling, and we expect Vi ~ If a/h^/m > 1 the instability acts on the neutrals but 
leaves the ions behind, and we expect large differences between Vi and Vn even at short 
wavelengths. 

These expectations are borne out by analysis. According to equation (21), the ampli- 
tudes of the ion and neutral displacements are A and B, respectively. Their ratio follows 
from equation (25), which for the problem at hand can be written in dimensionless form 
using equations (28) as 



The modulus of the relative amphtude is plotted as functions of wavelength in Figure 4 for 

the va > U and va < U cases, respectively. When the field is strong and h is large (short 
wavelength), the ions are held firmly in place and the instability is almost entirely in the 
neutrals, so is large. The numerical results are well fit by the asymptotic formula 

\B/A\ ~ /mh oc A^^. As h decreases the coupling improves, but the magnetic restoring 
force is still very large and the ion response never matches the neutrals. In this regime, 
\B/A\ ~ /mh^. Recall (see Figure 2) that in this limit the instability is of the weak form 
X ~ —iK^. 

When the field is weak, magnetic forces are relatively unimportant. In this case, \B/A\ 
is nearly unity over most of the range considered, although Figure 4 (bottom panel) shows 
a slight upturn at the shortest wavelengths, at which ion-neutral friction is relatively weak. 
The location of the upturn shows the expected scahng with B. For example, at A = 280, 
\B/A\ ~ 1 for S = l/zG, increasing to 25 if S = 10/iG and 50 if S = 2Q^lG. In all three 
cases, however, \B/A\ drops to unity for /i < 1 (A ~ 2.8 x 10^ AU). Not until the field is 



B_ 
A 



x'^ + h'^ — ix 



—IX 



(36) 
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strong enough to stabilize the bulk medium B ~ 230/xG does \B/A\ exceed unity for /i < 1. 

In deriving both the numerical and asymptotic results in the cases where x ~ —ih, the 
x'^ and /i^ terms in the denominator of eqn. (36) nearly cancel. Thus, the right hand side of 
the equation (36) must be computed using roots to the exact dispersion relation (equation 
(29)). 

The critical result from this analysis is that if S >~ lO/xG then at relevant length scales 
and timescales for bipolar outflows {I < 1000 AU and r < 1000 years), the fluctuations are 
primarily in the neutrals and not in the ions or magnetic field. In contrast to the MHD 
shear flow instabilities simulated by Jones et al. (1996) and Malagoli et al. (1996), the 
magnetic field cannot wind up and dissipate on small scales, and the instability is primarily 
hydrodynamic in character. We now turn to the observational implications of this conclusion. 

4. Observational Implications 

The calculations presented here do not permit quantitative predictions for observational 
signatures of the two instabilities for ambient ISM entrainmcnt by massive bipolar outflows. 
We would need to develop numerical models with two fluid MHD that incorporate non- 
planar geometry and non-linear growth, and would need chemical models and radiative 
transfer calculations to predict the resulting spectra. Such an analysis is beyond the scope 
of this paper. However, several physical signatures should be robust, and may be observable. 

4.1. Moleculair Spectroscopy 

As we showed in §3.3, unless B is less than about 10 /iG, over most of the parameter 
space of interest only neutrals participate in the resulting turbulence. That is, only neutrals 
will have a velocity significantly different from either the outflow jet or the ambient ISM. 
For example, if outflow jets had uniform velocity emerging from the protostar, then the 
emission from ionized molecules should not be detectable at velocities different from the 
outflow velocity. This signature should be particularly distinctive if there are ion species 
that are present in only one of these media. Outflow jets probably have a range in velocity, 
however, which will make isolating a velocity unique to the interface region difficult. 

If the transition between the jet and the ambient medium can be resolved, we might 
expect relatively large differences between the profiles of charged and neutral species. The 
ion flow would be laminar, and would show a smooth transition between the velocity of the 
core of the jet and the velocity of the ambient medium. The neutral flow would show the 



same range of velocities, but would be turbulent, leading to broader lines. 



4.2. Heating 

Figure 4 predict that throughout most of parameter space, the magnetic field prevents 
the ions from participating in unstable perturbations. The resulting velocity difference be- 
tween the ions and neutrals creates friction that can heat the gas. The heating rate per unit 
volume is given by 

T/ric = Pnl^ni{Vn " ^^^)^ ~ Pn^^m^^^ (37) 

For example, assuming I'n ~ 1 km s"^, n„ = 10^ cm~^, rii — 10~^ cm~^ gives F/^c ~ 6x 10~^^ 
ergs cm"^ s"^. Assuming that the heating is balanced by radiative coohng, the corresponding 
luminosity density is 0.04 Lq pc~^. The frictional heating rate dominates cosmic ray heating. 
Taking the latter from Spitzer (1978) we find 

= 1.5 X 10%vl (38) 

where is the neutral velocity in km/s. Unlike viscous heating, which is also strong in a 
turbulent boundary layer, F/^ic is independent of the thickness of the layer. 



4.3. Magnetic field configuration 

The orientation of the magnetic field in the outfiow could be mapped by measuring the 
polarization of infrared radiation emitted by the dust grains. According to our analysis, 
the magnetic field should remain well organized even though the neutral fiow is turbulent, 
provided that the mean field is aligned with the fiow (if the field were transverse to the fiow, 
the plasma would be unstable, and would carry the field with it, but the turbulent motions 
would be perpendicular to the field, and would not strongly bend it). 

Observations of a turbulent field in the outfiow could be explained if the ionization 
fraction were much larger than assumed here, resulting in stronger coupling, or if the grains 
were aligned by the turbulent neutral fiow rather than by the magnetic field. Observations 
of a straight field, however, demonstrate only that the field is well organized at or above 
the scale resolved by the observations (Heitsch et al. 2001). A low polarized intensity could 
be evidence for a disordered field below the resolved scale, but could also be due to poor 
alignment of the grains (Padoan et al. 2001). 
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5. Conclusions 

Outflows from massive protostars are observed to carry a large mass flux and to be 
poorly collimated. Both properties suggest that these outflows entrain substantial volumes 
of ambient ISM material. Turbulence driven by shear flow instabilities at the interface 
between the jet and the surrounding medium is one possible mechanism for entrainment. 

In this paper we have examined the Kelvin-Helmholtz instability in a weakly ionized 
medium with parameters chosen to be similar to molecular outflows. In order to study effects 
introduced by two fluids (plasma and neutral), we chose the simplest possible geometry: 
two uniform media (taken in all examples to have the same density) in relative motion 
at speed U, separated by a sharp interface, with a uniform magnetic field parallel to the 
flow. Under these conditions, a plasma is unstable if the relative velocity exceeds the Alfven 
speed (Chandrasekhar 1961). For the parameters of interest to our problem, the plasma is 
magnetically stabilized while the neutrals, in the absence of frictional coupling to the ions, 
would always be unstable. In the limit of perfect frictional coupling, the combined system - 
with the Alfven speed based on the total mass density instead of the plasma density - could 
be either stable or unstable, according to whether the bulk Alfven speed va exceeds U. 

Most of the analysis in our paper is based on the dispersion relation derived in §2, 
equation (30). Our main result is that for much of the size range of interest, the ions and 
neutrals are sufficently decoupled that the neutrals are unstable, while the ions are held in 
place by the magnetic field. At short scales, the growth rate is well approximated by the 
growth rate of the hydrodynamic instability, kU. At longer scales, the fastest growing mode is 
either the magnetically modified MHD mode, if the magnetic field is weak enough (va < U), 
or, if va> U , new mode with growth rate kU {kU /v). The transition between these modes 
of instability occurs at kU/v ^ 1. The dispersion diagram is plotted in Figure 2. The 
relative amplitudes of the ion and neutral fluctuations are shown in Figure 4, demonstrating 
that over most of the regime of interest, the ions drop out of the instability. 

We beheve these results continue to hold in cylindrical geometry. In fact, they should be 
even more pronounced because at short wavelengths the cylindrical and planar cases coincide, 
while at long wavelengths the instability is suppressed by geometry. We also expect that 
more general perturbations, such as those with a large transverse wavenumber ky, would be 
characterized by similar ion-neutral decoupling. 

Although we have not followed the instability into the nonlinear stage, there is no reason 
why frictional coupling should be more efficient at large amplitude, so we expect our results 
to apply to the nonlinear regime. There is, however, a compelling case to be made for 
nonlinear calculations: they are necessary to produce detailed observational predictions. As 
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such, they should follow the thermal and chemical state of the gas, as well as its dynamics. 

Based on the linear theory, however, it appears that our model may be tested in several 
ways, mentioned in §4; differences in line profiles between charged and neutral species, 
efficient heating by ion-neutral friction, and the polarimetric signature of a well organized 
magnetic field. We will report on spectroscopic evidence in a forthcoming paper (Watson et 
al. 2004). 

Our results do, however, leave us with a conundrum. Churchwell (1997) pointed out 
that the rate of entrainment required to explain the mass flux in outflows is unreasonably 
large, if the entrainment is hydrodynamic in nature. The low degree of ion-neutral coupling 
predicted by our model suggests that the entrainment is more hydrodynamic than mag- 
netohydrodynamic, although nonlinear simulations and analysis are necessary to settle the 
issue. 

We are happy to acknowledge support by NSF grant AST-0328821, NSF PHY-0215581 
and by the Graduate School of the University of Wisconsin, Madison. F.H. is grateful for 
support by a Feodor-Lynen Fellowship from the Alexander von Humboldt Foundation. 

A. Ions Held Stationciry 

We begin with the same assumptions as specified in section 2, with the additional 
assumption that the ions remain at rest. That is, they are not disturbed by the perturbation. 
The equation of motion for the neutral particles is: 

Pndt^n + PnV„ • Vv^ = "VPn " Pn^niyn (Al) 

Assuming a periodic perturbation, we obtain: 

Pn{0J + Ukx)Un = -kxPn + iPnl^Un (A2) 
Pn{uJ + Ukx)Wn = id^Pn + ipnl^Wn- (A3) 

Taking the curl and combining to eliminate the pressure term, we obtain 

Pn{u; + Ukx-ti^){dl-kl)w^ = (A4) 
We integrate over the interface to obtain, 

A{pn{uj + Uk)d,Wn} = iA{pni^d,Wn} (A5) 
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We introduce a solution of the form 



Wni = A{u^Uk)e^' z<0 (A6) 

Wn2 = A{uj + Uk)e-^' z>0 (A7) 

to obtain 

Pn{oJ + Ukf -ipni^iuj -Uk) = (A8) 
After converting to dimensionless notation, we obtain 

x'^ + h'^-ix = (A9) 

^(i± V-1 - 4/i2) (AlO) 

Ui±i{l + 2h')) (All) 



X 



X 



For h<l, 

X f« i, -i/i^ (A12) 

Thus, we see that the —ih'^ mode observed in section 2 originates from the uniform drag 
that the neutral particles feel from the uniformily distributed ionized particles. 
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Table 1. Unstable Roots of Characteristic Equation 



x"^ — imx^ — a^x"^ — i{mh'^ - 


- a'')x - a^h^ = 


(1) (2) (3) (4) 


(5) (6) 


Root Dominant Terms 


Physical Regime 


—ih 3, 6 


h>l 


-ih"^ 5, 6 


1> h and > h 


i{\±h) 3,4,5,6 


-^>h>l 
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Fig. 1. — Diagram of the standard Kelvin-Helmholtz instability. Since the fluid is assumed 
to be incompressible, a small periodic perturbation at a shear layer will tend to grow. In 
general, the vorticity perturbation {STy) is out of phase with the vertical displacement of the 
vortex sheet (rj) . If it can be arranged that the vorticity has the distribution shown by the 
dashed line, then the equilibrium flow sweeps negative vorticity toward point P. The effect 
is as if an "extra eddy" , rotating counterclockwise were added to the flow. The centrifugal 
force associated with the eddy pushes the interface between P and Q outward, resulting in 
the growth of rj with time. 
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Fig. 2. — The growing instability root of the characteristic equation as a function of wave- 
length. The strong magnetic field regime (top) assumes U =10 km s^^ and B =lmG and 
the weak magnetic field regime (bottom) assumes U =10 km s~^ and B =l/xG. The crosses 
represents the numerically determined root. The solid lines represent the analytic approx- 
imation to root, h (h) at short wavelengths and h {h?) at long wavelengths for the strong 
(weak) magnetic field regime. 



Growth Rate in Different Rcqir-ios 
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Fig. 3. — The growth rate of the Kelvin-Helmholtz instabihty depends on the density, flow 
velocity and magnetic field strength. The transition between the standard Kelvin-Helmholtz 
instability (growth rate oc h) and the new, slower instability (growth rate oc h^) is given for 
3 magnetic field strengths. 
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Fig. 4. — The amplitude of the neutral perturbations relative to the ion perturbations. Nu- 
merical results are represented by crosses. The lines represents the analytic approximations 
to the relative perturbations as specified above. The strong magnetic field regime (top) as- 
sumes U —10 km s~^ and B — 1 mG and the weak magnetic field regime (bottom) assumes 
U =10 km s-^ and S = 1 /xG. 
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